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ABSTRACT 


Integration of the equations of rotation of the planet Mercury 
is facilated by having a computer perform the pertinent numerical 
algorithms . The perturbations of the two thirds orbital period 
rotation of Mercury are calculated accurate to one tenth of an arc 
second. The numerical solutions to the libratural motion of Mercury 
are presented. 



In a recent report 


3 


we presented the analysis of the resonance- 


locked rotation of the planet Mercury under the influence of the 
planetary potential and the orbital eccentricity. The purpose of the 
present note is to record the numerical solution of the librational 
motion. Integration of the equations of rotation of the planet Mercury 
is facilated by having a computer perform the pertinent numerical 
algorithms. The perturbations to the two thirds orbital period 
rotation of Mercury are calculated accurate to one tenth of an arc 
second. 

The equation of the planetary rotation about its polar axis 
has the form 

= f-r 


where f is the true anomaly and CP, denoting the angle between the 
smallest of the moment of inertia with the radius vector, is the 
solution of the following differential equation^ 





A 



in which measures the difference between the two 

small momentsof inertia and e is the orbital eccentricity. 

By applying the Runge-Kutta integration algorithm, equation ( 2 ) 

is solved on an IBM 7090 computer. The computer program performs the 
necessary algorithms for successively updating expansions of 
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equation (2) truncating only insignificant terms. The over-all pre- 
cision is consistent with the 10 7 truncation level. 

In order to investigate the instantaneous physical librations of 
the planet Mercury, the variable of function $ has been transformed 
from the true anomaly f to the mean anomaly M. If E is the planet's 
eccentric anomaly, a well-known formula of orbital motion 



( 5 ) 


constitutes the basis of the transformation. Then, by Kepler's equa- 
tion the mean anomaly can be calculated according to 




Solutions of § are generated for various combinations of para- 
meter X and initial conditions for a period of 100 years . For the 
orbital eccentricity of Mercury, the value e = 0.206 is adopted. 
The results of the period and the small amplitude of this libra- 
tional motion are: 


x = 0.00000$ x = 0.000050 x = 0.000500 

Period Amplitude Period Amplitude Period Amplitude 
(years) (days) (years) (days) ( years ) ( days ) 

80 0.008 25 0.008 8 0.008 
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For large amplitude libration, the period should be corrected by 
a factor of an elliptic integral. The results indicate that it wquld 
be possible to determine the parameter X for the planet Mercury if the 
period and the amplitude of its physical libration can be observed. 

Since the resonance condition is due to the angle jr 
about — -~r Jtf — C> , the rotation of the planet Mercury can remain 
locked-in motion when the rotational rates are within the range be- 
tween an upper and a lower limit. A series of computer runs estab- 
lished the following results: 


X 

Lower limit for 
rotational period 

(days) 

Upper limit for 
rotational period 

(days) 

0.000005 

58.538 

58.75^ 

0.000050 

58.260 

59.032 

0.000500 

57.381 

59.911 
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